Derealization of relativistic Dirac particles in disordered one-dimensional systems 

and its implementation with cold atoms 
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We study theoretically the localization of relativistic particles in disordered one-dimensional 
chains. It is found that the relativistic particles tend to derealization in comparison with the 
non-relativistic particles with the same disorder strength. More intriguingly, we reveal that the 
massless Dirac particles are entirely delocalized for any energy due to the inherent chiral symmetry, 
leading to a well-known result that particles are always localized in one-dimensional systems for 
arbitrary weak disorders to break down. Furthermore, we propose a feasible scheme to detect the 
derealization feature of the Dirac particles with cold atoms in a light-induced gauge field. 
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Since the pioneering work by Anderson[l|, substantial 
efforts have been devoted to the understanding of trans- 
port properties of electrons in disordered systems Q. A 
very significant advance along this direction is the scaling 
theory proposed by Thouless et al [3| and Abrahams et 
al 0| ■ In the scaling theory, it is argued that the quantity 
j3 = din g/dhiL is a monotonic and nonsingular function 
of g only, with g a dimensionless conductance and L the 
sample size. Currently there exists a well-known result 
that the conductance g approaches to zero as the sam- 
ple size L goes to infinite for any disordered one dimen- 
sional systems. In particular, an arbitrary weak disorder 
strength leads to the localization of all states of electrons 
in one-dimensional chains [1, [H, 0] . 

Actually, an implicit pre-condition for the above re- 
sults is that the particles are governed by the Schrodinger 
equation since the (quasi-)particles addressed in con- 
densed matter systems are in general non-relativistic. 
Notably, relativistic Dirac particles have recently been 
paid a significant attention because the quasi-particles in 
honeycomb lattices (such as electrons in the graphenepj 



and cold atoms in the optical lattices 



H|) .ultra-cold 



atoms in a light-induced gauge field [11|, [12, 
and trapped ions [l|| may be described by the relativis- 
tic Dirac equation. However, Anderson localization in the 
relativistic region has been less studied in literature 17 1. 
Apart from the Klein's finding[l8| that the transmission 
of Dirac particles is essentially different from that of non- 
relativistic particles, it is also fundamentally important 
and interesting to study the aforementioned localization 
issue of relativistic particles and to work out how to sim- 
ulate the predicted results with currently available tech- 
niques. 

In this paper we study Anderson localization of rela- 
tivistic particles in disordered one-dimension systems by 
using the finite scaling method and the transfer-matrix 
technique. The finite-size scaling analysis reveals that all 
the states of the massive relativistic Dirac particles are 



localized in the systems; while, the localization length of 
relativistic particles is longer than that of non-relativistic 
particles with the same disorder strength. More intrigu- 
ingly, the states of the massless Dirac particles are en- 
tirely delocalized for arbitrary disorder strength due to 
the chiral symmetry, providing a distinct example that 
breaks down the well-known result for non-relativistic 
particles that all states are localized in disordered one- 
dimensional systems. Furthermore, we propose how to 
simulate the wanted Dirac particles in disordered one- 
dimensional chains and the related (de)localization prop- 
erties with recently /newly developed techniques in the 
cold atomic systems 0, HH, E IH HI ■ 

Let us consider a particle transmitting through a one- 
dimensional chain with N rectangular barriers as shown 
in Fig. 1(a), where the potential 



V(x) = 







others 



(n = l,2,--.,A0 (1) 



with V n being a constant randomly distributed in the 
range [—5,5]. Here S represents the disorder strength. 
For simplicity, we assume that x' n — x n = a and x n +\ — 
x' n = d for any n. A relativistic particle with the mass 
m and energy E is inserted into the N barriers from the 
left, which is described by the Dirac equation 

\~ihca x 4- + mc 2 a z + V(x) - EU(x) = 0, (2) 
ax 

where c denotes the velocity of light, a x , z are the Pauli 
matrices, and ip(x) represents a two-component spinor. 
A general solution of Eq. ([2]) for any region with a con- 
stant potential, e.g., V n in Fig. 1(a), is given by 



4>{x)=A[ I )e^+B[ 1 )r ,:-!) 



where p n represents the momentum of the particle, n n — 
(E - mc 2 - V n )/(cp n ) and (E - V n ) 2 = m 2 c 4 + c 2 p 2 . If 
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E and V n are fixed then p n can in principle be either 
positive or negative. Here we choose the positive one, 
and thus the coefficients A and B denote the amplitudes 
of the spinors moving along the positive x-axis and its 
opposite direction, respectively. 

We now look into the transmission for N potentials 
as shown in Fig. 1(a). Denoting the amplitudes of the 
spinor at a position approaching to x n (x' n ) with the 
infinitesmall amount from the left (right) as {A n ,B n } 
({A' n ,B' n }), we may obtain a relation between the am- 
plitudes based on the continuity of the wave-function, 
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where denotes the transfer matrix of the n-th bar- 
rier, and its elements are given by 



(AOn^cos^ 



(M, 
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(5) 



c 2 p 2 + m 2 c 4 . 



with k = [E — mc 2 )/(cp) and E 2 

For comparison, we recall the results for non- 
relativistic case which is described by the Schrodinger 
equation [-^-£2 + V(x) - E k ]$(x) = 0, where E k is 
non-relativistic kinetic energy. In this case, we have a 
relation similar to Eq. but M„ is replaced by M„ 
which is the transfer matrix for the nth barrier calcu- 
lated by the Schrodinger equation. It is straightforward 
to derive that the elements of the matrix M„ have the 
same form of Eq. (|5]), but K, p n (and K n ) are replaced 
by the non-relativistic counterparts p s = \Z2mEk and 
pf t = y/2m(E k - V n ), respectively. 

We consider an incoming relativistic particle with en- 
ergy E (Ek for non-relativistic kinetic energy ) from the 
left, then the amplitude of the outgoing at the right side 
of the N barriers is related to the amplitude of the in- 

A r \ , r J I A[ 



coming by the relation ^ ^ j — M I ^ j , where the 

total transfer matrix M J (hereafter we use the super- 
scripts J = D and S to denote the relativistic and non- 
relativistic cases, respectively) for the N barriers reads 



M J = M J N D J M 



L N-l 



■ D J MlD J Mf. 



(6) 



Here D J = diag{exp(— ip J d/K), exp(ip J d/h)} represents 
the displacement matrix between two nearest neighbor 
barriers, with p D — -J E 2 /c 2 — m 2 c 2 . 

The transport properties can be extracted from the 
transfer matrices for both relativistic and non-relativistic 
cases. At zero temperature, the conductance through 
the N barriers is given by the Landauer formula [2 5j 
G J = ^.g J , where g J = l/\(M J )u\ 2 is the dimen- 
sionless conductance. It is noted that the localization 
length £ J or the Lyapunov exponent -f J is defined as 



7 J = 1/£ J = — liniL^oo (In g J )/L, where L is the total 
length of the chain L = Nb = N ( we choose b = d + a 
as the unit of length), and (• • •) denotes the averaging 
over the disorders. The £ J is a function of the energy E 
(E r ) and can be used to characterize a localized state: a 
state is a localized state if £ J is finite and is a delocalized 
(extended) state if £ J is divergent. 

It is hard to obtain an analytical expression for the 
localization length £ J in a general case; however, it has 
been shown that limL^ 00 ((ln5' s ) / L) always exists for any 
energy of a non-relativistic particle in an arbitrary weak 
one-dimensional disordered system In a similar way, 
we can find that limi^ oc ((ln g D )/L) also exists for rel- 
ativistic massive particles. The numerical procedure for 
both relativistic and non-relativistic cases is as follows: 
one can define a J N = W ln 9i( N ), where g-(N) 

is the conductance for a specific configuration of fixed 
N barriers, and then a J N is an averaged quantity for the 
number of N c configurations. We find that, for both non- 
relativistic and relativistic massive particles, a conver- 
gent value a J N can always be derived for sufficient large 
N, which could be considered as the localization length 
£ J of the system. The result implies that the state for 
the massive Dirac particles is also a localized state for 
arbitrarily weak disorders, as in the non-relativistic case. 
The localization length (Lyapunov exponent) as a func- 
tion of the potential width a is plotted in Fig. 2. It is seen 
that the localization length of the relativistic particles is 
always longer than that of the non-relativistic particles. 

We now turn to examine the validity of the single- 
parameter scaling equation (3 J = in the relativistic 
region. The scaling quantities (3 J as a function of (In g J ) 
for m — 2.5 x 10~ 4 mo with mo the mass of a nuclear in 
both relativistic and non-relativistic cases are plotted in 
Fig. 2(b). The non-relativistic case was studied in Ref.0], 
and our results are essentially the same as those presented 
there. It is seen clearly from Fig. 2(b) that the assumption 
of the single-parameter scaling theory is still valid for the 
massive Dirac particles in the disordered systems. 

It is interesting to note that an analytical result for 
a massless particle (m = 0) can be derived, from which 
one is able to find unexpectedly that the particle is en- 
tirely delocalized for arbitrary energy. For massless parti- 
cles, the transfer matrix of total N barriers is obtained as 
M = diagje^/' 1 , e~^/ h } with ip = -Npb + J2n =1 p n a. 
The transmission amplitude t — exp(i<p/h) is a pure 
phase factor, and the dimensionless conductance g D = 1. 
In this case, the localization length £ D for the massless 
particles approaches to infinity, and thus breaks down the 
famous conclusion that the particles are always localized 
for any weak disorder in one-dimensional systems. 

The inherent physics is simply the chiral symmetry. 
The time-independent Hamiltonian for the Dirac parti- 
cles is given by Hp = —ihca x -^+mc 2 a z + V(x), and the 
chiral operator for a Dirac spinor is the matrix 7 s = a x in 
one-dimensional cases. Under the discrete chiral trans- 
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formation the spinor is transformed as ip c 
transformed Hamiltonian 



H c = Y°H D 1 5 = -ihccr x - mc 2 a z + V(x). (7) 

ax 

Then the chirality is conserved for a massless particle. 
Noting that j 5 <fi± — ±</>± with cf>± = f _^ , the gen- 
eral solution of the massless Dirac equation described in 
Eq.© can be rewritten as ij}{ x ) = A(j) + e^ px +B(j)-e~^ px , 
i.e., the first (second) term is actually the eigenstate of 
the chiral operator with positive (negative) chirality. As- 
suming that the incoming wave function is ipi n {x) = 
A<f> + efi px , then the outgoing wave-function ip ou t( x ) — 
A'(f> + e% px + B'4>-e~% px should have the same positive 
chirality, i.e., the reflection rate B' must be zero for mass- 
less particles because of the conservation of the chirality. 
However, the chirality is not conserved for a massive par- 
ticle, so the reflection B' in principle can not be always 
zero. In this case the massive particle should be localized 
for any weak disorders. Alternatively, an intuitive picture 
of localization of relativistic particles may be understood 
with the Klein's paradox. For massive particles, when 
the height of the potential barrier reaches the order of 
mc 2 , the barrier becomes nearly transparent. Since some 
of the barriers are transparent, the localization length in- 
creases. For the massless particles mc 2 is zero and thus 
every barrier is transparent. Therefore, massless parti- 
cles are always delocalized. 

We now turn to address how to simulate the rela- 
tivistic particles with cold atoms. We consider the adi- 
abatic motion of atoms having a tripod level configu- 
ration in the field of three laser beams, as shown in 
Fig.l(b,c)[H[Il[Iil. The ground states |1),|2) and |3) 
are coupled to an excited state |0) through spatially vary- 
ing laser fields, with the corresponding Rabi frequencies 
Qi,Q,2 an d ^3 j respectively. The full quantum state of the 
atoms $(r) can be described as $(r) = Ylj=o 0j ( r ) |j) , 
where r is the atomic position. The original Hamiltonian 

p 2 

of the atom with the mass m a takes the form H = — h 
V H (r) + V L (r) + H int , where V H (r) = £- =0 ^(r)|j)(j| 
represents an external harmonic trapping potential, and 
Vt(r) denotes a state- independent random potential. 
Hi n t is the laser-atom interaction Hamiltonian given by 
Hint = — ?l 23 J= i (fij j 0) | + H.c.), where the Rabi fre- 
quencies are chosen as fti = ft sin 9e~ %kx /^/2, ft2 = 
ft sin 6>e' fc: 7\/2, and ft 3 = ft costfe"^. Here ft = 
\/|fti| 2 + |ft2| 2 + |ft3| 2 , k is the laser wave vector and 
the angle 9 defines the relative intensity 13, 13 , 3] • Fol- 
lowing Ref. 12(, one is able to obtain an effective one- 
dimensional Dirac-type Hamiltonian as 



7 5 V> and the sin 4 9, c^ — cos 9 is the effective 'speed of light'. 

In the derivation, we have assumed that the trapping 
potential V H (r) is independent on the internal states. 
Comparing the original Dirac Eq.([2]) with the Dirac- like 
Eq.© achieved in cold atoms, the effective 'speed of 
light' in cold atoms is c* and the effective mass m — 
^ tan 2 #sin 2 9. Note that the mass m of the simulated 
Dirac particle is not the mass m a of the cold atom itself 
and it is a remarkable feature that the mass m in the 
simulated Dirac-likc equation can be controlled by the 
laser beams. Thus both the massive and massless Dirac 



H k 



C+CTxPx + lz<Jz + V"(x) + V L {x), 



(8) 



equations can be realized with cold atoms [27 1. 

Finally, we discuss briefly the detection of the local- 
ization length. For concreteness, we assume that the 
tripod-level configuration in Fig. 1(b) is provided by the 
atoms of 87 Rb, where the ground states (Jl), |2), |3)} 



are the hyperfine states 5 ^/2^ = l,mp 



-1,0,1) 



up to an irrelevant constant, provided that the wave 
vector of the atoms p x /h <C fccos0f26|. Here ^ z 



and the excited state |0) is given by either the state of 
5 2 P 3 / 2 (F = 0) or 5 2 P 3/2 (F = 2,m F = 0). In this case, 
the harmonic potential Vff(r) has been experimentally 
realized by the far-off resonant laser beams in the imple- 
mentation of the spinor condensates of 87 Rb 24( . In addi- 
tion, a feasible approach to detect the localization length 
of the relativistic particles can follow the case in non- 
relativistic particles implemented in Ref. [2(| , except that 
three additional laser beams represented by ftj are re- 
quired. The experiment starts with an elongated cluster 
of ultracold 87 Rb atoms trapped by a harmonic potential 
Vh- A far-off-resonance laser beam (such as wavelength 
1.06/im used in Ref.(2fj|) creates an optical wave guide 
along the horizonal x axis, and a loose longitudinal trap 
is also realized by such laser beam. Three laser beams 
with the resonant wavelength of rubidium ( wavelength 
0.78 /im, near the resonant also works) but different po- 
larizations, as shown in Fig. 1(c), are shined on the atoms 
to create an atomic gas that could be described by the 
Dirac equation. The longitudinal confinement is switched 
off at time t — 0, and the atomic gas starts to expand 
in the guide along the x direction. A disordered poten- 
tial Vl is applied to the expanding atomic gas using an 
optical speckle field produced by passing a laser beam 
(with the wavelength about 0.514 /iim[2(j) through a dif- 
fusing plate. One then detects the spatial distribution of 
the atoms at increasing evolution time using absorption 
imaging. As in the experiments 1|J 2^| , one can directly 
measure the localization length of the particles by the 
density profiles. Therefore, the comparison of Anderson 
localization between relativistic and non-relativistic cases 
can be made for the conditions with and without the ad- 
ditional laser beams ftj. Considering that the two main 
ingredients to observe (de)localization of Dirac particles, 
the disordered potentials [2(| and light-induced gauge 
field [23)], have been achieved in the recent experiments 
on the atoms of 87 Rb, it is expected that the cold atoms 
may offer a novel platform for the study of Anderson lo- 
calization in the relativistic region. 
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In summary, we have found that the relativistic parti- 
cles tend to delocalization and revealed that the mass- 
less ones are entirely delocalized in disordered one- 
dimensional systems. The predicted features may be 
tested by future experiments with ultra-cold atoms. On 
the other hand, the (de)localization of the relativistic 
particles may also be observed in a disordered graphenc, 
where Dirac electrons are confined to move in one dimen- 
sion and the impurities are small enough such that the 
scattering does not occur between the two Dirac points. 

We thank L.-M. Duan for many fruitful discussions. 
This work was supported by the RGC of Hong Kong, 
the NSFC (No. 10674049), and the SKPBRC (Nos. 
2006CB921800 and 2007CB925204). 



potential, which is sufficiently strong along the transverse 
direction of the x axis, such that the original massive 
atoms are actually, or at least in principle, confined in 
a one-dimensional guide along the x axis. A crucial con- 
dition of this Dirac equation is p a /h ~ l/l a <C kcosO 
(a = x,p), where l x (l p ) is the longitudinal (transverse) 
confinement and k ~ 8.0 x 10 6 m -1 for the atoms of 87 Rb. 
If we simply take 9 = arccos(l/8), this condition implies 
that the dynamics of the atoms along the transverse di- 
rection is nonrelativistic for l p < 1.0 fim and thus the 
atoms may still be confined by the original trap, while it 
is relativistic along the x axis for l x S> 1.0/xm. 
[27] If we choose a specific potential given by Vi H — V3 1 = 
h 2 k 2 sin 2 8/2m a , we are able to obtain a massless one- 
dimensional Dirac equation as the 7 Z in Eq. vanishes. 
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FIG. 1: Schematic representation of the system, (a) N rect- 
angular potentials, (b) Atom with tripod level structure in- 
teracting with three laser beams, (c) The configuration of the 
laser beams to realize a Dirac-like equation with the lasers Qj 
and a disordered potential with the laser Ql- The atoms are 
confined in a one dimensional wave guide by a harmonic trap. 




FIG. 2: (a) The localization length £ J as a function of the 
potential width a. The insert is the corresponding Lyapunov 
exponent 7 J .(b) The scaling quantity /3 D (/3 s in the insert) as 
a function of (\ng D ). (\ng J ) is the average of 1000 configura- 
tions. The other parameters are E 
with the units of mc 2 . 
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